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CRITICAL SHEAR STRESS OF AH 


LONG ELATE IN THE PLASTIC REGION 


By Elbridge Z. Stowell 


SUMMARY 


The plasticity reduction factor in the formula for the buckling 
stress of a long plate under uniform shear has been computed from the 
properties of the stress-strain curve for the material. Some limited 
tests on the shear buckling of 2kS-0 aluminumr-allpy plates in the 
plastic region tend to confirm the value of the reduction factor 
predicted by the theory . Restraint against rotation along the edgeB 
of the plate has little or no effect upon the value of the reduction 
factor according to theoretical calculation. 


INTRODUCTION 


The problem of computing the critical shear BtreBS for an infinitely 
long plate with edges elastically restrained aga i nst rotation was solved 
for the elastic stress range in reference 1. The theory of reference 1 
is no longer applicable , however, in cases where buckling occurs beyond 
the elastic range; in such cases, the theory must be modified to allow for 
curvature of the stress— strain curve. The problem consists in the compur- 
tation of the plasticity reduction factor rj in terms of the properties 
of the stress-strain curve for the material. 

The general equations for the buckling of a thin plate under combined 
compression and shear in the plane of the plate were derived in reference 2. 
Upon the assumption that the plate remains in the plastic state during 
buckling, the differential equation of equilibrium of the plate was derived, 
together with the correspo nding expressions for the variations in energy at 
buckling. These general expressions were then applied to the case of a 
plate in which one of the compressive forces and the shear force were zero, 
that is, to a plate under simple compression. 

The same general expressions may also be applied to the case of a 
plate under pure shear in the plane of the plate by setting both camr- 
pressive forces eq ual to zero. This application has been made in the 
present paper. 
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RESULTS 


Theoretical results.— The plate is assumed to he infinitely long 
•with side edges elastically restrained against rotation. The applied 
shear stress is assumed to he uniformly distributed and , as in 
reference 2, no unloading of the plate is assumed to occur during 
buckling . By virtue of these assumptions, the plate remains in the 
plastic state during buckling and the relations developed in reference 2 
appby. 

The procedure followed is to employ the same deflection surface in 
the energy expressions for the plastic case that was used in reference 1 
for the elastic case. From the resulting expressions the approximate 
value of the critical shear stress in the plastic region is computed. 

The expression for the critical shear stress in such a case is 


Jt 2 ksEt 2 

*12 Cl - |i 2 )b 2 


in which 

r critical shear stress 

k s coefficient depending upon conditions of edge restraint 

and shape of plate 

t thickness of plate 

b width of plate 

E Young’s modulus 

M- Poisson’ s ratio 

t] coefficient which allows for reduction in shear stress in 

plastic range 

The value of the wave length required to make the critical shear stress 
a minimum is found and is substituted in the expression for the critical 
shear stress, which is then known in terms of the angle of the waves, the 
magnitude of the elastic restraint, and the plasticity coefficients. 
Division of thin expression for the critical shear stress by the corre- 
sponding approximate expression in the elastic range gives the value of 
the coefficient rj to a closer approximation than either of the indi- 
vidual critical stresses. 

The coefficient t} is given formally in the appendix as a function 
of the angle of the waves the magnitude of the elastic restraint e. 



fl-nrt the plasticity coefficients E a /E and C^. The resulting equation is 


n = ¥ 


E 0 sto 20 o 2 \/l 2^ + 2[l + 2 sin20 - (l - 03) cos^Jf^e) 


2\/fl(0 +2(1 + 2 sin%)f 2 (e) 


where 


1 _l 

— 00 a x 


secant modulus of material 


fl(e), f 2 (e) functions of elastic restraint e 
C3 plasticity coefficient 

0O wave angle in elastic range 


(3 ~ C3)f 2 (e) 

~2\/f l(«j / 


1 - C, 


jin^20 


+ (3 + 03)f 2 («) 


The angle of the waves 0 is that angle -which vill make the critical shear stress a 

minimum. As soon as the magnitude of the elastic restraint is selected, the angle of 

the waves, and therefore tj, “sy computed for any stress. 

Computations of tj were made for a plate of 2 ljS-T aluminum alloy in the oases -where 
the edges of the plate were simply supported and clamped. The procedures used are given 

in detail in the appendix. The calculations with both edge conditions resulted in 

E„ E + 

substantially the same curve, shewn in figure 1 . Curves of tj - and r| ■=■ — + as 

3E 35 

obtained for compression in reference 2 are included for c carp arisen* 
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Experimental check .— The critical shear stress of eight long 
plates of 2ljS— 0 al uml mim alloy has "been measured in the plastic 
region hy Gerard (reference 3) • An attempt vas made to secure a 
clamped condition along the edges end buckling wan determined "by the 
difference in str a i n on opposite sides of the sheet as measured with 
electrical gages. Gerard’s data are represented hy the test points 
in figure 2. The solid— line curve of figure 2 gives the theoretical 
valueB of i) as computed hy equation (l) from the axial stress- 
strain curve supplied hy Gerard for 24S-0 al uminu m alloy in reference 3. 
This theoretical curve confirms the trend of the experimental points in 
a satisfactory manner, although it lies in the upper part of the scatter 
hand. 


In reference 3 Gerard suggested that the shear secant modulus he 
used as the effective modulus for a long plate under shear. Gerard’s 
assumption is, from equation (2) of reference 3, 



where the G’s refer to the slope of a shear stress— strain curve 
derived from the axial stress— strain curve. The dash— line curve in 
figure 2, which was computed from figure 8 of reference 3, ■ shows the 
value of t] that results from this suggestion. The agreement with 
the solid— line curve, which represents the values of t) computed from 
the theory of the present paper, is seen to he fair. 

The necessary conditions for use of the shear secant modulus, 
however, should he examined further. In the appendix, equation (l) 
is shown to reduce to the form 

T] = Constant 

if C3 is constant. Eor Gerard’s curve for 24S— 0 alumlmrm alloy the 
requirement for a constant C3 is closely fulfilled. When the value 
of p is computed from equation (l), it is found that 

i - °- 8 9 ?r 

This relation is plotted as the solid line in figure 2. Besults of 
the calculation are also given in a more familiar form in figure 3. 

Caution should he employed in generalizing the use of the shear 
secant modulus to materials, such as 2te3-T al umin um alloy, for which 
C^ is not a constant in the plastic range. The upper curve in 
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figure 1 represents E g /E for 2te-T aluminum alloy and therefore 
describes the axial stress— strain curve for that material. The center 
curve, which represents T|, cannot he obtained from the upper curve by 
multiplying the ordinates by a constant factor throughout the plastic 
range; therefore r, for this material is not the shear secant modulus. 


CONCLUSIONS 


The theory of plastic buckling previously reported has been applied 
to the determination of the plasticity reduction factor in the formula 
for the buckling stress of a long plate under uniform shear from the 
properties of the stress— strain curve for the material. Some limited 
data on the shear buckling of 2kS-0 aluminum-alloy plates in the plastic 
region tend to confirm the value of the reduction factor predicted by 
the theory. The theory indicates that restraint against rotation along 
the edges of the plate has little or no effect on the value of the 
reduction factor. 


Langley Aeronautical Laboratory 

Hational Advisory Committee for Aeronautics 
Langley Field, Ya., May 2k , 19^3 
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APPENDIX 

ANAIXSIS 


General theory.— Reference 2 shows that for a plate under shear 
alone the work T done on the plate, eapressed in a Cartesian (z*, y') 
system, is 


T = Tt [l 


where 


T 

applied shear stress 

w 

deflection of plate at (z 


half-wave length 


The expression for the strain energy in the plate Y-^, from 
reference 2, is 




[ 8^w ^ 2 8^w 

Ww Bz « 2 v 2 



(3) 


where 



1 5t 

2E S 



E t 


tangent modulus of material 


If, in addition, the plate has equal elastic restraints e against 
rotation along the side edges, the strain energy in these restraints Y 2 
is, from reference 2, 
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Y 2 = 


D‘e 

2b 


1 A. 

2 r ^ 12 


.1 

2 2 


dx* + 


D* € 
2b 


X 

2 

.JL 

2 


cm 


L^rTJy. = . 


± 

2 


dx* 


(b) 


Let a more suitable oblique coordinate system be defined "by 


x = x* + y* tan 0 


y ’ 

^ ~ COB <jS 


(5) 


as In reference 1. (See fig. 4.) 

Then T, Yi, and Y 2 i* 1 new system are 


T = Tt 


h, p L 

2 j 2 

LI I i_k 


dw dw /dw\ 2 
dx dy + 


(ID Bln ^ 


dx dy 


( 6 ) 


Yt = 


D* 


L, 

2 / 2 


2 008 ^ j_£l J_A 
2 2 


cos 




- 2 (I - C 3 ) Bin 2 ^ 



dx 2 y 


_i ¥ (p) 2 + [ 1 + 4t^_ (1 _o 3) ](^) 


+ [l + 2 tan 2 ? -( 1 - 03)] 32 " S2,r 


+ tan 0 d^w 
cos (f dxdy 


dx 2 dy 


^ + & - (1 - Co) cos 2 ^ 

dx 2 dy 2 3 dx 2 


J 


dx dy 


(7) 
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V 2 = 


D»e 

2b cos 2 ^ 



^ + I 31 sin 
ay 02 




+ 


bw 

bx 



(8) 


where 

. _ h 

1 cos $ 

A comparison with the corresponding formulas in appendix B of reference 1 
will show that the two groups are identical if Poisson* s ratio is taken 
as 0.5 and if C^ - 1; that is, if only the elastic range is considered. 

As in reference 1 , the assumed deflection -surface is 


cob 


in which the half-^wave length A will be adjusted to make the critical 
stress a minimum. Substituting this deflection surface into 
eapjressiops (6), (7), and (8) and performing the indicated integrations 
gives 



T 



jt^bpt sin 0 
2 A 


Lt 


Vi = w Q 2 2^^ — 

1 0 Ub-jA cos $ ] 


f— 1 2 (1 - C^ sin^ll + M 

' [cos 2 ^ ' • J / hi cos 0' 

\ A 


1 


+ 2 1 + 2 Bin g ^ - (1 - C 3 ) cos 2 ^ ^ 


cos 




( 9 ) 


To = w. 


2 rt 2 D*eA 
2bi^ cos30 
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where 


M = (s - ^) £2 + (5 - ^) 6 + 1 


Setting up the equality T = Y-j_ + Y 2 , which holds at buckling, and 
solving for the critical shear stress t gives 


r = 


sin 2$ 


r i E 

— “W — 2 ( 1 “ c q) sin 2 ^ + — — ~t • p 

j^cos 2 ^ ' J /t>x COB 0^? 


2 + 2 a ^- D -^V — (l — C 3 } cos 2 ^ n 


2e 


cos 




cos f h-^^t 


or, replacing cos 0 "by "b, and writing 
and ^ = f2(^) yields 


M + 


2c 

* 2 


= fx(e) 


Sta ^] + 7^ 

w 

+ 2 [jL + 2 sin 2 ^ — (1 — C-^ cos 2 ^ f 2 (e) 




rt 2 !)* 

h 2 t sin 2$ 


( 10 ) 
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In the elastic range (C^ = l) this expression reduces to -the corresponding 
formula in reference 1. 

The half-wave length X and the wave angle 0 will adjust themselves 
to make t a minimum; that is. 


8t _ 

which gives 

/hf \fa.(Q 

^ \Jl Sin 2 ^ 


cos 






from, which 


1 _ 1 


min 


sin 20 


2 ^Te) \jl - i-^3. sin^20 


+ 2 |l + 2 sin 2 ^ — (l — cos r 2 («) 


l ah’! 

h^t 


( 12 ) 


The angle 0 is found from the relation 


<dr 


which gives 


cos 20 = 


= 0 

00 


(3 ~ 03 ) f 2 (g) 


2\/f^U) 

, — ■ ■ - - + (3 + Co^ f£(e) 


(13) 



The value of 0 la thus fixed as soon as the restraint is chosen and the value of E^/!E B is 

assumed. This value of 0 la the one to he used In equation (12) for computation of the 
minimum critical shear stress 


In the elastic range, the corresponding minimum critical shear stress (to) , 1b obtained 

TTn n 

by putting Cr3 = 1 and D’ = D = in equation (12): 

o 9 


(tq) 


JL 


min sin 20 q L 


2\(fi(e) + 2 (l + 2 sin%) f 2 (e ) 


'-1 b 2 t 


(ib) 


■where the subscript 0 signifies that the quantity is in the elastic range. Thus, from 
equation (13)., 


cos 20 o 


fg(e) 


\Jf p(e) + 2f 2 ( € ) 


( 15 ) 


The value of q 
stress (tq) 

ml u 


is the ratio of the critical shear stress r ^-n to the oritloal shear 
that would he obtained if the material were wholly elastic; thus 


T] 


E fl sin 20 o 2 V 


(l - — ■ 2 °-^ sin 2 ^ + 2 

[l + 2 sin 2 ^ 

— | 

(l “ 0 3 ) cos 2 0]f 2 (O 

® Bin 20 


2 Vfl(«) + s( 

L + 2 sin 2 0q) 

If2(e ) 


The value of q as a function of stress 'will now he considered for the two extreme 
oases of simply supported and olamped edges: 
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(a) For the case of simply supported edges, 

f x (e) = f 2 ( € ) = 1 

cos 2^0 = 3 
sin 20 o = 0.943 

sin% = 0.333 
cos 20 = 


and 


3 ~ c 3 


\jl - 23. sin 220 


+ 3 + C 3 


E s 0.943 2 \Jl — - — g C 3 aln.220 + 2 [l + 2 Bin 2 0 - (1 - C 3 ) cos 2 0] 

5^33 


r, = 

1 E sin 20 


This value of T) is plotted against stress in figure 1 for 24S-T aluminum 
alloy. 

(t>) For the cane of clamped edges, 

6 =« 

f x (e) = 5-1^ 
f 2 (e) = 1.24 
cos 2 0 q = 0.26 
sin 20o = O.966 
sin^Q = 0.37 


cos 20 = 


1 » 2 ^(3 ~ C3 ) 


4.53 


]jl — — g C 3 sin g 20 


+ 1.24(3 + C3) 
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and 


T| = IS- 


_ 0.966 
1 sin 20 


4.53 \j 1 5-^ sin^ + 2.48 [l + 2 sin 2 ^ - (1 - C 3 ) cos 2 ^] 


8.84 


This value of tj is nearly identical numerically -with the value from 
case (a) and is represented by the same curve in figure 1. 

Condit ions under which Gerard’s shear secant modulus applies.— 
Equation (1) gives the general formula for -which involves, 

besides E b /E, the following quantities: 


EH 

Functions only of restraint 

BEE 

coefficient e 

0o 

Angle of waves on assumption of 
elasticity (by equation (15)) 
depends only on 6 

0 

Angle of waves in plastic range 
(by equation (13)) depends on e 
and plasticity coefficient C3 

c 3 

Plasticity coefficient 


The magnitude of the restraint coefficient e is considered to he 
part of the data furnished with the problem. Thus, fq( e ) and f2( e ) 
are known constants; the angle 0q is therefore also a known constant. 
If, in addition to these quantities entering into equation (l), the 
plasticity coefficient C^ should also be a constant, equation (l) 
would reduce to the form 


E 

T| = Constant r~ 


Gerard’s assumption is, from equation (2) of reference 3, 
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where the G* s refer to the slope of a shear stress-strain curve derived 
from the axial stress— strain curve. 


The plasticity coefficient C^ is defined as 
the condition for constancy of C 3 is that 


1 + 1 

2 + 2 Eg 


and therefore 


E t =XE S 


or 



x£i 

e i 


where X is same constant. The relation = Ke-j*, where K is a 

constant, satisfies this requirement. Eor Gerard* s curve for 2kS—0 aluminum 
alloy, X is approximately 0.28 over the plastic range, and the requirement 
for a constant Cq is closely fulfilled. When the value of q iB computed 

from equation (l) it is found that q = O .89 -^r. This relation is plotted 

as the solid line in figures 2 and 3 . 
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Figure I.- Values of 77 for a long plate of 
24S-T aluminum alloy under uniform 
shear stress. (The curves for E s /E and 
Efy/E for compressive stress c r aire 
included for comparison.) 
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Figure 2.- Values of 77 for a long plate of 
24S-0 aluminum alloy under uniform 
shear stress. 
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